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We study the magneto-elastic coupling behavior of paramagnetic chains in soft polymer gels exposed to 
external magnetic fields. To this end, a laser scanning confocal microscope is used to observe the morphology of 
the paramagnetic chains together with the deformation field of the surrounding gel network. The paramagnetic 
chains in soft polymer gels show rich morphological shape changes under oblique magnetic fields, in particular a 
pronounced buckling deformation. The details of the resulting morphological shapes depend on the length of the 
chain, the strength of the external magnetic field, and the modulus of the gel. Based on the observation that the 
magnetic chains are strongly coupled to the surrounding polymer network, a simplified model is developed to 
describe their buckling behavior. A coarse-grained molecular dynamics simulation model featuring an increased 
matrix stiffness on the surfaces of the particles leads to morphologies in agreement with the experimentally 
observed buckling effects. 


I. INTRODUCTION 

Magneto-responsive hybrid gels (MRGs) have been attract¬ 
ing great attention due to their tunable elasticity, swelling 
properties and shape that can be remotely controlled by a mag¬ 
netic field. They have potential applications as soft actuators, 
artificial muscles, as well as sensors OHa and can serve as 
model systems to study the heat transfer in hyperthermal can¬ 
cer treatment 0. Compared to other stimuli-responsive gels, 
MRGs have the advantage of fast response, controlled me¬ 
chanical properties and reversible deformabilities Had. A 
typical MRG consists of a chemically cross-linked polymer 
network, swollen in a good solvent, and embedded magnetic 
particles 0111 . The size of the magnetic particles can range 
from ~ 10 nm to several jim Q. 

The origin of the magneto-responsive behavior of MRGs 
is the magnetic interaction between the magnetic filler parti¬ 
cles as well as their interaction with external magnetic fields 
mmol. In a uniform magnetic field, paramagnetic particles 
can be polarized and act as approximate magnetic dipoles. 
Depending on their mutual azimuthal configuration, the dipo¬ 
lar interactions can be either attractive or repulsive. In a liq¬ 
uid carrier, the dipolar interaction drives the magnetic parti¬ 
cles to form chains and columns MM aligning in the di¬ 
rection of the magnetic field. However, in a polymer gel, the 
particles cannot change their position freely. Instead, rela¬ 
tive displacements of the particles, induced e.g. by changes in 
the magnetic interactions, lead to opposing deformations of 
the polymer network. As a result, the magnetic interactions 
can induce changes in the modulus of the gel QUa. This 
magneto-elastic effect is well known to be related to the spa¬ 
tial distribution of the magnetic particles MM- For exam¬ 
ple, the modulus of anisotropic materials that contain aligned 
chain-like aggregates of magnetic filler particles ifTSi l22l - E4l 
can be significantly enhanced when an external magnetic field 
is applied along the chain direction ||2l. The anisotropic ar¬ 
rangement of particles also dominates the anisotropic magne¬ 
tostriction effects MM- 


Different theoretical routes have been pursued to investigate 
the magneto-elastic effects of MRGs; macroscopic continuum 
mechanics approaches l29l . mesoscopic modeling 1113 - 
fT9l . and more microscopic approaches Il30l - l32l that resolve 
individual polymer chains. Theoretical routes to connect and 
unify these different levels of description have recently been 
proposed ll3TH35l . The authors of Ref. |M] show how the in¬ 
terplay between the mesoscopic particle distribution and the 
macroscopic shape of the sample affects the magneto-elastic 
effect. In addition to these factors, recent experiments and 
computer simulations also point out that a direct coupling be¬ 
tween the magnetic particles and attached polymer chains can 
play another important role 0 13011311 [3313911 . 

An experimental model system showing a well-defined par¬ 
ticle distribution and a measurable magneto-elastic effect can 
help to understand the magneto-elastic behavior of MRGs 
at different length scales. Projected into a two-dimensional 
plane, the distribution of magnetic particles in thin diluted 
MRGs can be detected using optical microscopy or light scat¬ 
tering methods ESI Sol. By combining these techniques with 
magnetic or mechanical devices, it is possible to observe the 
particle rearrangement when the MRG sample is exposed to 
a magnetic field or mechanical stimuli ESI Ell- For three- 
dimensional (3D) characterization. X-ray micro-tomography 
has been used li23l . Here we introduce another 3D imag¬ 
ing technique - laser scanning confocal microscopy (LSCM). 
Compared to normal optical microscopy, LSCM is able to ob¬ 
serve 3D structures and it has a better resolution Ea. Com¬ 
pared to X-ray micro-tomography, LSCM is faster in obtain¬ 
ing a 3D image and easier to combine with other techniques 
for real-time investigation IS3E1- 

We use LSCM to study the magneto-elastic effects of para¬ 
magnetic chains in soft gels. As a result, we find that the 
paramagnetic chains in soft gels (elastic modulus < 2 Pa) un¬ 
der an oblique magnetic field show rich morphologies. De¬ 
pending on the length of the chain, modulus of the gel and 
strength of an external magnetic field, the chains can rotate, 
bend and buckle. The deformation field in the polymer net¬ 
work around the deformed paramagnetic chains can also be 
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tracked. The result confirms that the chains are strongly cou¬ 
pled to the polymer network. A simplified model is developed 
to understand the magnetically induced buckling behavior of 
the paramagnetic chains in soft gels. In addition to serving as 
a model experimental system for studying the magneto-elastic 
effect of MRGs, our approach may also provide a new mi- 
crorheological technique to probe the mechanical property of 
a soft gel na. Furthermore, our results may be interesting 
to biological scientists who study how magnetosome chains 
interact with the surrounding cytoskeletal network in magne¬ 
totactic bacteria Il46l . 


II. MATERIALS AND METHODS 

The elastic network was obtained by hydrosilation of a 
difunctional vinyl-terminated polydimethylsiloxane (vinyl- 
terminated PDMS, DMS-V25, Gelest Inc.) prepolymer 
with a SiH-containing cross-linker (PDMS, HMS-151, Gelest 
Inc.). Platinum(O)-1,3 -di vinyl-1,1,3,3-tetramethyIdisiloxane 
complex (Alfa Aesar) was used as a catalyst. A low- 
molecular-weight trimethylsiloxy terminated PDMS (770 
g/mol, Alfa Aesar GmbH & Co. KG, in the following “PDMS 
770”) served as a solvent that carried the polymer network and 
the paramagnetic particles. The paramagnetic particles were 
purchased from microParticles GmbH. They were labeled 
with fluorophores (visible in LSCM). The materials consist 
of porous polystyrene spheres. Within the pores, nanopar¬ 
ticulate iron oxide was distributed rendering the particles su- 
perparamagnetic. To prevent iron oxide leaching, the parti¬ 
cles had a polymeric sealing that also held the fluorophores. 
The particles had a diameter of 1.48 ±0.13 /rm HTl . We 
measured the magnetization curve izi of the paramagnetic 
particles by a vibrating sample magnetometer (VSM, Lake 
Shore 7407). We found about 20% deviations in the mag¬ 
netic properties of the magnetic particles (e.g., magnetic mo¬ 
ment Il47l '). In order to observe the deformation field in the 
polymer network, we used fluorescently labeled silica parti¬ 
cles as tracers. They had a diameter of 830 ±50 nm and the 
surface was modified with N,N-dimethyl-N-octadel-3-amino- 
propyltrimethoxysilylchloride. 

The paramagnetic particles were dried in a vacuum oven at 
room temperature overnight before they were dispersed into 
PDMS 770. The prepolymer mixture was prepared with 9.1 
wt% vinyl-terminated PDMS and 90.9 wt% SiH-containing 
cross-linker. The prepolymer mixture (2.86 wt%) was dis¬ 
solved in PDMS 770, which contained the paramagnetic par¬ 
ticles. Finally, by adding PDMS 770, which carried the cat¬ 
alyst, the concentration of the prepolymer mixture in the sol 
solution was carefully adjusted in the range from 2.74 wt% to 
2.78 wt%. This concentration range guaranteed the formation 
of soft gels with an elastic modulus lower than 10 Pa. In the 
sol solution, the catalyst concentration was 0.17 wt%, and the 
concentration of magnetic particles was 0.09 wt%. The sol 
solution was agitated at 2500 r/min with a Reax Control (Hei- 
dolph, Schwabach, Germany) for 2 min for homogenization, 
followed by ultrasonication (2 min. Transsonic 460/H, Elma) 
to disperse the magnetic particles. Then the sol solution was 


introduced into a thin sample cell (~ 160 /rm thick and ^ 1 cm 
wide) by capillary forces. The sample cells consisted of two 
No. 1 standard coverslips, separated by ^ 160 /rm spacers. 
After sealing with two-component glue, the cells that con¬ 
tained the sol were exposed to a 100.8 ±0.5 mT magnetic 
field. The paramagnetic particles aligned into chains along 
the direction of the applied magnetic field while the prepoly¬ 
mer was crosslinking. A visible reaction of the prepolymer 
occurred within 10 min, and the rheological measurements 
showed that it took about 40 min to form a gel. Due to the 
low concentration of magnetic particles, the magnetic chains 
in the gel were well separated (> 30 /rm). The length of the 
chains varied from a single particle up to about 170 particles. 
In some samples, 3 wt% silica particles were added as trac¬ 
ers. We stored the samples at ambient temperature for at least 
two weeks before testing. The modulus of the gel in the sam¬ 
ple cells was characterized using microrheological techniques 

Ell mi. 

A home-built LSCM setup was used to observe the chain 
structure in the gel Il43ll4^ . We were able to analyze a sample 
of thickness of about 150 fim. A homogeneous magnetic field 
was attained by building Halbach magnetic arrays near the 
sample stage of the LSCM ||49]| . A 32-magnet array (Fig.^) 
was used to change the field direction while keeping the field 
strength constant (216.4 ± 1.1 mT BtI I. Another 4-magnet 
Halbach array l47l was used to change the field strength (up 
to 100.8 ± 0.5 mT). The magnetic field was measured by a 
Lake Shore Model 425 Gaussmeter with a transverse probe. 


III. RESULTS 

In the absence of a magnetic field, the paramagnetic chains 
in a soft gel kept the aligned morphologies Il47l . When a mag¬ 
netic field (216.4 ± 1.1 mT) was applied in the direction par¬ 
allel to the chains (Fig. [^, images for 0°), the paramagnetic 
chains still aligned with the original chain direction (horizon¬ 
tal). We changed the direction of the magnetic field step-by- 
step (5°) in the clockwise direction (^1 min between steps, 
quasi-static). We also define the orientation of the magnetic 
field B as the angle included between the magnetic field and 
the initial chain direction (see Fig. [^). The left images of 
Fig# show a short chain with 15 particles in a gel of stor¬ 
age modulus G' of 0.25 ± 0.06 Pa. The chain rotated to fol¬ 
low the magnetic field. However, the rotation angle of the 
chain is smaller than the orientation angle of the magnetic 
field (Fig.[g3). This difference increased until the orientation 
of B reaches 135°, where the chain flipped backward and had 
a negative angle. The chain again became parallel to the field 
when the orientation ofB increased to 180°. The morphology 
of the chain was the same at orientations of the magnetic field 
of 0° and 180° because of the superparamagnetic nature of the 
particles. Note that the chain was not straight at the interme¬ 
diate angles (e.g., images for 60°, 90° and 120°). Instead it 
bended. 

The images on the right-hand side of Fig.[T}; show a longer 
chain with 59 particles in the same gel. When the orientation 
of B was 30°, the chain rotated and bended, with its two ends 
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FIG. 1: (a) Laser scanning confocal microscopy (LSCM) was used 
to observe the morphologies of the paramagnetic chains in the soft 
gels. The Halbach magnetic array provided a homogeneous mag¬ 
netic field (here B = 216.4± 1.1 mT). This array could be rotated to 
change the orientation of the magnetic field, (b) The orientation of 
the magnetic field B was successively increased from 0° to 180° in 
36 steps (square points). A magnetic chain of 15 particles rotated to 
follow the magnetic field, but the rotation angle was smaller than the 
orientation angle of B (dashed line), (c) Morphologies of magnetic 
chains in a soft gel change when the orientation angle of B increased. 
The scale bar is 10 /Im. The gel in (b) and (c) had a storage modulus 
G' of 0.25 ±0.06 Pa. 


tending to point in the direction of the magnetic field. When 
the orientation of B was 60°, the chain wrinkled and started 
to buckle. A sinusoidal-shape buckling morphology was ob¬ 
served when the magnetic field was perpendicular to the orig¬ 
inal chain (orientation of the magnetic field of 90°). When the 
orientation of B increased from 90° to 120°, the left part of 
the chain flipped downward in order to follow the magnetic 
held. The right part flipped upward when the orientation of B 
increased from 120° to 150°. Finally, when the held direction 
was again parallel to the original chain direction (orientation 
of the magnetic held of 180°), the chain became straight. The 
same rotation/buckling morphologies as in Fig. could be 
observed when increasing the orientation of B from 180° to 
360°. 

We also directly applied a perpendicular magnetic held to 
the paramagnetic chains in the soft gels. The paramagnetic 
chains showed different buckling morphologies (Fig. |^) de¬ 
pending on the chain length. Fig.|^ gives frequency counts of 
the different morphologies in the same sample (G' = 0.25 ± 
0.06 Pa) under a magnetic field of 100.8 ± 0.5 mT. In total 
180 chains were observed. Longer chains tended to buckle 
with a higher number of half waves. Moreover, the distribu¬ 
tions overlapped, implying that the paramagnetic chains with 
the same length could have different morphologies under the 
perpendicular magnetic held. 

These buckling morphologies are reminiscent of the buck¬ 
ling of paramagnetic chains in a liquid medium under a per- 



FIG. 2: (a) Different morphologies of paramagnetic chains in a 
soft gel {G' = 0.25 ± 0.06 Pa) under a perpendicular magnetic field 
(100.8 ±0.5 mT). The original chain direction was horizontal, and 
the applied magnetic field was vertical. The scale bar is 10 pm. 
(b) Frequency count of different buckling morphologies in the same 
sample. M is the number of half waves. 


pendicular magnetic field ISOllSTl . The most stable morphol¬ 
ogy in the latter system was a straight chain aligning along the 
magnetic field direction. However, in our system this mor¬ 
phology was not observed. Even the short chains showed 
a rotation angle smaller than the orientation of the magnetic 
field (e.g.. Fig. [^). The major difference between our exper¬ 
iments and Refs. l50l and ISTI was the nature of the surround¬ 
ing medium. In our system, the polymer network around the 
paramagnetic chains impeded the rotation of the chains into 
the magnetic field direction (a more detailed discussion will 
be given below). 

We used IMAGE! software (NIH 1521) to extract the skele¬ 
tons of the chains which have 2 half waves (S-shaped). The 
amplitude of deflection or deformation of different chains was 
quantified by the square root of the mean square displacement, 
i.e. Amplitude = ((y^) — where y measures the par¬ 

ticle displacement along the field direction. The results are 
shown in Eig. The amplitude increased with increasing 
chain length. At the same chain length, the amplitude tends 
to increase with increasing magnetic field strength (Eig. [^; 
an example is also given in Eig. |^) or with decreasing gel 
modulus (Eig.[^). 

The modulus dependence of the amplitude demonstrated 
that the polymer network around the paramagnetic chains im¬ 
peded the chain deformations. Therefore, the deformation 
field within the polymer network plays an important role to 
understand the buckling of the chains. We thus added tracer 
particles into the gel matrix, and used their trajectories to 
record the deformation field around the paramagnetic chains. 
As shown in Eig. a linear paramagnetic chain buckled and 
formed an S shape in a perpendicular magnetic field. The am¬ 
plitude increased with increasing field strength, while the con¬ 
tour length of the chain remained constant. The chain exten¬ 
sion decreased along the original chain direction (horizontal 
direction) and increased along the perpendicular direction. Si¬ 
multaneously, the polymer network around the chain followed 
























4 



Chain length (gm) 



FIG. 3: Influence of chain length, strength of magnetic field and 
elastic modulus of the gel matrix on the amplitude of the S-shaped 
chains, observed when the magnetic field is applied perpendicu¬ 
larly to the initial chain orientation, (a) The elastic modulus of 
the gel was 0.25 ± 0.06 Pa, and the magnetic field strengths were 
216.4± 1.1 mT (hlack squares), 80.5±0.4mT (red triangles) and 
18.7 ±0.1 mT (blue circles), respectively, (b) The magnetic field 
strength was 216.4 ± 1.1 mT and the elastic moduli of the gel were 
0.015 ± 0.005 Pa (black squares), 0.25 ± 0.06 Pa (red triangles) and 
0.78 ± 0.22 Pa (blue circles), respectively. The solid lines represent 
linear fits and are included as guides to the eye. 


the deformation (Fig. [^) of the paramagnetic chain, both in 
the transverse and longitudinal directions. This confirmed that 
the paramagnetic chain is strongly coupled to the polymer net¬ 
work. Within our experimental resolution, the chain seemed 
to have a rigid non-slip contact to the surrounding network. 


IV. MODELING 



FIG. 4: (a) Influence of the magnetic field strength on a buckling 
chain. From top to bottom, the magnetic field strengths were 0 mT, 

10.1 ± 0.1 mT, 26.4± 0.1 mT, 38.6 ± 0.2 mT and 60.2 ±0.3 mT, re¬ 
spectively. The modulus of this gel was about 0.01 Pa. (b) Tracer par¬ 
ticles were inserted into the gel matrix of the sample. Tracking these 
embedded tracer particles, the deformation field in the gel matrix was 
determined (the Particle Tracker plug-in developed on IMAGE! soft¬ 
ware was used for this purpose) (54|. The red solid line represents 
the skeleton of the magnetic chain shown in (a) for a field strength of 

60.2 ±0.3 mT, and the dashed blue line indicates the original chain 
shape, (c) We modeled the paramagnetic chain in the elastic gel as 
a continuous object uniformly carrying dipolar magnetic moments. 
Without the magnetic field, the straight chain was oriented along the 
x-axis. Under a perpendicular magnetic field B (oriented along the 
y-axis), the magnetic chain deformed. The surrounding polymer net¬ 
work impeded the chain deformation. 


We now turn to a qualitative description of the situation in 
the framework of a reduced minimal model. Theoretically 
capturing in its full breadth the problem of displacing rigid 
magnetic inclusions in an elastic matrix is a task of high com¬ 
plexity and enormous computational effort 15^ . We do not 
pursue this route in the following. Instead, we reduce our 
characterization to a phenomenological description in terms 
of the shape of the magnetic chain only. This is possible if the 
dominant modes of deformation of the surrounding matrix are 
reflected by the deformational modes of the magnetic chain. 

Below, we assume N identical particles on the chain. In the 
undeformed state, the straight chain is located on the x-axis of 
our coordinate frame. The contour line of the deformed chain 
running through the particle centers is parameterized as y(x), 
see Fig.[^. 


A. Magnetic Energy 


integral and shifting the path of integration from the contour 
line of the chain to the x-axis, we obtain the magnetic interac¬ 
tion energy im 

m 1 

Emagn=W I . dx, (1) 

where xi and X 2 label the end points of the chain. The prefac¬ 
tor W has the dimension of energy per unit length and is given 

bylEI 


W 


47td^ ’ 


( 2 ) 


where /Tq is the vacuum magnetic permeability, m the mag¬ 
netic moment of a single particle, and d its diameter. 


First, concerning the magnetic energy along the chain, we 
assume dipolar magnetic interactions between the particles. 
In the perpendicular geometry, the external magnetic field ap¬ 
proximately aligns all dipoles along the y-axis. For simplicity, 
we only include nearest-neighbor magnetic interactions. In an 
infinite straight chain, this would result in an error given by 
a factor of i^(3) « 1.2, where ^ is the Riemann zeta function 
Il33l|55]|56l. Within our qualitative approach this represents 
a tolerable error. Replacing the magnetic interaction energy 
between the discrete magnetic particles by a continuous line 


B. Elastic Bending Energy 

Next, we need to include terms that provide a measure for 
the magnitude of the elastic deformation energy. To estimate 
the importance of different modes of the elastic matrix defor¬ 
mation, we analyze the experimentally determined displace¬ 
ment field around the distorted chain shown in Fig. |^. For 
this purpose, we model the continuous matrix by a discretized 
spring network IISlEil. Network nodes are set at the posi¬ 
tions where the displacement field was tracked experimentally 
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FIG. 5: The four most occupied normal modes of the deformation 
in Fig. & after projection to an elastic spring network, ordered by 
decreasing magnitude of contribution to the overall deformation. 
The normal modes (a), (b), and (d) are of an “oscillatory” type, 
whereas mode (c) represents a longitudinal contraction. Correspond¬ 
ing relative weights of the modes are = 0.095, = 0.057, 

= 0.055, and = 0.051, where we normalized the sum of the 
weights over all modes to unity. For better visualization, the over¬ 
all amplitudes are rescaled as against the actual weights. The matrix 
region in close vicinity of the chain is indicated by black arrows. 


by tracer particles. The nodes are then connected by elastic 
springs. After that, we determine the normal modes of de¬ 
formation of this network Et). Finally, we can decompose 
the experimentally observed deformation field in Fig. & into 
these normal modes. Occupation numbers give the contri¬ 
bution of the «th mode to the overall deformation. 

The four most occupied modes are shown in Fig.|^ We find 
a major contribution of “oscillatory” modes, i.e. alternating 
up and down displacements along the central horizontal axis. 
Such oscillatory displacements of the matrix are induced or at 
least reflected by corresponding oscillatory displacements of 
the chain, see Fig.|^. A bending term of the form BTI 

Ebend = Cb r ( 3 ) 

[l+y(;c)2]5/2 

becomes nonzero when such deformational modes occur and 
is therefore taken as a measure for their energetic contribu¬ 
tion. In addition to that, we have experimental evidence that 
the chain itself shows a certain amount of bending rigidity 
in, possibly due to the adsorption of polymer chains on the 
surfaces of the magnetic particles. Similar indication follows 
from two-dimensional model simulations, see below. 


C. Elastic Displacement Energy 


therefore any displacement of an inclusion induces a distor¬ 
tion of the surrounding gel matrix. We model this effect by a 
contribution 


Edispi = Q r [y{x)f [1 +yix)^] dx. (4) 

Jxi 

This term increasingly disfavors the rotations of longer 
straight chains, which reflects the experimental observations 

El. 

Moreover, in Fig. the third dominating mode of the ma¬ 
trix deformation corresponds to a contraction along the chain 
direction and an expansion perpendicular to it. We conjec¬ 
ture that this should be the dominating mode in the deforma¬ 
tional far-fleld, yet this hypothesis needs further investigation. 
It is induced by chain deflections in y-direction, which imply a 
shrinking extension in .r-direction (experimentally we observe 
that the chain length is conserved under deformations and that 
the individual magnetic particles remain in close contact). We 
simultaneously use Eji^pi to represent the energetic contribu¬ 
tion of this type of underlying matrix deformation. 


D. Energetic evaluation 

We now consider the resulting phenomenological model 
energy E,ot = E^agn + Ebend + Edispi- First, we only address 
the bulk terms of the energetic expressions. Minimizing them 
with respect to the functional form of y{x) and linearizing the 
final expression, one obtains, above a certain threshold of the 
magnetic held, wave-like oscillatory deformations BtII . This 
is in agreement with the observation of the wrinkles at onset 
in Fig. [TJ; and the final oscillatory shapes in the inner part of 
the longer chains in Fig. [^. 

Detailed knowledge about the boundary conditions of the 
deflection and its derivatives at the end points of the chain 
would be necessary to fully determine the chain shape from 
the above equations. These boundary conditions depend on 
the interaction with the matrix and are not accessible in the 
present reduced framework. Therefore, we proceed in a dif¬ 
ferent way. From the experiments, the shape of the chains 
is known and can to good approximation be represented by a 
polynomial form 


M-l 

y{x) = s n {x — mb) ioxx\<x<X 2 ^ (5) 

m=0 

where M is again the number of half-waves, the prefactor 
S sets the strength or amount of chain deformation and de¬ 
flection, b is the spacing between the nodes, and the inter¬ 
val [xi,X 2 ] follows from the experimental result of preserved 
chain length L, 


J yJl + \y'ix)ydx = L. 


( 6 ) 


The bending term does not energetically penalize rotations 
of a straight chain, see Fig. for M = 0. Yet, such rotations 
cost energy. Boundaries of the block of material are fixed. 


The polynomial form reproduces the experimentally observed 
straight chain ends as well as the smaller oscillation ampli¬ 
tudes inside longer chains. 
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FIG. 6: Contributions to the total energy as a function of the strength 
5 of deformation and minimized with respect to xi and X 2 for a chain 
of the shape given by Eq. jg. Here we show the cases (a) M = 2, 
L = 3b and (b) M = 4, L = 4.5b. The total energy E,ot has a global 
minimum as a function of S, which corresponds to the most stable 
chain shape. We always observed the global minimum for symmetric 
shapes. 


Next, we insert Eq. Q into the above expressions for the 
energy and minimize with respect to 5, jci, and X 2 for a given 
M, with the constraint of constant length L, see Eq. (j^. Pa¬ 
rameter values of the coefficients Cb and Cd are found by 
matching the resulting shapes to the corresponding experi¬ 
mental profiles (chain deformations for G' = 0.25 Pa and mag¬ 
netic field B = 100.8 mT as in Eig.l^, M = 2, are used for this 
purpose). We obtain Ct « O.OlWlr and Q « 2W jb^. 

To illustrate how the energetic contributions vary under in¬ 
creasing preset deformation, we plot in Eig.j^the energies for 
increasing S for two fixed combinations of M and L. The total 
energy E,ot shows a global minimum in both panels, which 
we always observed for symmetric chain deformations. As 
expected, with increasing amplitude S the magnetic energy 
decreases, whereas the deformation energies increase. 

Next, we determine the minimal total energy as a function 
of chain length L for increasing number of half-oscillations M. 
We can see in Pig. that with increasing chain length L the 
shapes that minimize the energy show an increasing number 
of half-waves M in good agreement with the experimental data 
in Pig. 1^. 

Moreover, we quantify the amplitude of the chain deflection 
or deformation by 

/ ^ ^ • dx 

Amplitude = J (/) - (jf ^ (.) = —!-. (7) 

'' X2-X1 


FIG. 7: Energies of chain deformations of the shape given by 
Eq.fg, minimized with respect to 5, xi, and X 2 , as a function of chain 
length L and number of half-oscillations M. Each curve describes a 
shape of M half-waves with a minimum total length of (M — 1 )fo. The 
resulting curves show crossing points from where the total energy for 
an increasing L is lowered by bending one extra time (jumping to a 
higher M) rather than conserving the same shape. 

Resulting values are plotted in Fig. and compared with cor¬ 
responding experimental data. As mentioned above, we op¬ 
timized the model parameters with respect to the experimen¬ 
tal data for a magnetic field intensity of B = 100.8 mT. We 
demonstrate in Fig.j^that moderate variations of the magnetic 
field intensity only slightly affect our results: the brighter 
curves are obtained when multiplying the magnetic energy 
scale W by a factor ~ 1.42, cotTesponding to an increased 
magnetic field intensity of approximately B ^ 216 mT BtI . 
This is in agreement with the experimental observations. We 
include in Fig. the experimentally determined values for 
B = 80.5 mT and B = 216.4 mT. Only a slight trend of in¬ 
creasing deflection amplitudes is found for this increase of 
magnetic field intensity. 

Together, although the curves for M = 2 in Fig. slightly 
overshoot the data points, Figs.|7]and[^are in good agreement 
with the experimental results. The amplitude of deflection 
and deformation is not observed to unboundedly increase with 
chain length L in the experiments. Fikewise, our model pre¬ 
dicts that longer chains prefer to bend one extra time (switch¬ 
ing to higher-M shape) rather than to show too large deflection 
amplitudes. 


V. COARSE-GRAINED MOLECULAR DYNAMICS 
SIMULATION 

We also studied the buckling of the chain using two- 
dimensional coarse-grained molecular dynamics simulations 
by means of the ESPResSo software Il58ll59l . A simple model 
was developed that allowed us to analyze the influence of par¬ 
ticular interactions and material properties on the buckling ef¬ 
fect. Here, we focus on the elasticity of the polymer matrix in 
the immediate vicinity of the magnetic particles. 

By choosing the coarse-grained scale for our model, we 
ignore any chemical details but rather describe the system 
in terms of the magnetic particles as well as small pieces 
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FIG. 8: Resulting deflection amplitudes of the chain deformation, 
calculated according to Eq. 0- Darker curves represent the model 
parameters optimized with respect to the experimental shapes for 
a magnetic field intensity B = 100.8 mT. Brighter curves were 
obtained by increasing the magnetic energy scale IT by a factor 
~ 1.42, which corresponds to an increased magnetic field intensity of 
g~216mTpi7l. comparable with the triangular experimental data 
points. Both, model curves and included experimental data points, 
demonstrate that moderate variations of the magnetic field intensi¬ 
ties only slightly affect the observed deflection and deformation am¬ 
plitudes. The value of b necessary to perform the analysis was de¬ 
termined from the M = 2 experimental data as b = 12.6 /tm. For 
M >2 “bumps” appear in the curves, which arise from a change in 
the type of deformation as illustrated by the insets: for each M curve, 
the end points of the chain for lower L correspond to nodes of the 
deformational oscillations, i.e. y(x\) ~ y{x 2 ) — 0 (lower left inset); 
for higher L, these outer nodes shift to the inside of the chain (upper 
right inset). As seen from Fig.|^ the lower-L parts of the curves are 
not energetically preferred. 


of polymer gel. As the buckling effect appears to be two- 
dimensional, and as the ground states for systems of dipolar 
particles have also been found to be two-dimensional ll5^ . 
we use this dimensionality for our simulations. We study a 
chain of 100 magnetic particles with a significant amount of 
surrounding elastic matrix. 

As in the analytical approach, the gel matrix is modeled 
by a network of springs. Here, however, we use a regular 
hexagonal mesh as a basis. To mimic the non-linear elastic 
behavior of polymers, we use a finitely extensible non-linear 
elastic spring potential (FENE-potential ||60|) for the springs 
along the edges of the mesh. As a simple implementation of 
the finite compressibility, we introduce EENE-like angular po¬ 
tentials with a divergence at zero and 180 degrees on the an¬ 
gles at the mesh points BtII . The magnetic particles are mod¬ 
eled as rigid spheres interacting by a truncated, purely repul¬ 
sive Lennard-Jones potential, the so-called Weeks-Chandler- 
Andersen potential MM- Their magnetic moment is as¬ 
sumed to be determined purely by the external magnetic field 
and to be constant throughout the simulation. I.e. we assume 
that the external field is significantly stronger than the field 
created by the particles. The magnetic moments are taken 
parallel to the external field and with a magnitude given by 
the experimentally observed magnetization curve. The cou¬ 
pling between the particles and the mesh is introduced in such 
a way, that under the volume occupied by a particle, the mesh 
does not deform, but rigidly follows the translational and ro¬ 


tational motion of the particle ll47l . Hence, a local shear strain 
on the matrix can result in the rotation of a magnetic particle, 
but not its magnetic moment. 


An important point is the elasticity of the polymer matrix in 
the immediate vicinity of the magnetic particles and, in partic¬ 
ular, between two magnetic particles. We study two situations 
here, the first one including a stiffer region in the immediate 
vicinity of the particles, the second one without such a stiffer 
layer and directly jumping to the bulk elasticity. The stiffer 
layer, if imposed, is created using a spring constant larger 
by four orders of magnitude on those springs which originate 
from mesh sites within the particle volumes El. The angular 
potentials are unchanged. 


A comparison between the cases with and without a stiffer 
layer of gel around the magnetic particles can be seen in 
Fig. 0 The images show a small part of the resulting config¬ 
uration of magnetic particles and the surrounding mesh for a 
field applied perpendicular to the initial chain direction. Thus 
the magnetic moments of the particles are oriented perpen¬ 
dicular to the undistorted chain direction. This results in an 
energetically unfavorable parallel side-by-side configuration 
for the dipole moments. The energy can be reduced either 
by increasing the distance between the dipoles along the ini¬ 
tial chain direction, or by moving dipoles perpendicularly to 
the initial chain direction so that they approach the energeti¬ 
cally most favorable head-to-tail configuration. If the matrix 
is made stiffer immediately around the particles, and thus the 
contour length of the chain cannot change significantly, the re¬ 
positioning towards the head-to-tail configuration causes the 
buckling effect seen in the experiments (Eig. |^. When one 
assumes the matrix immediately around the magnetic parti¬ 
cles to be as soft as in the bulk of the material, neighboring 
particles can move apart and the chain breaks up into individ¬ 
ual particles or small columns perpendicular to the original 
chain direction. Additionally, a layer of increased stiffness 
also introduces a bending rigidity of the chain. In Eig. 10 the 
full chain and the surrounding matrix is shown for an exter¬ 
nal field of magnitude 216 mT, which from the experimental 
magnetization measurements corresponds to a magnetic mo¬ 
ment of about 4.5 • 10^*^ Am^ ll47ll . Due to the different di¬ 
mensionalities, the elastic modulus of the surrounding matrix 
could not be directly matched to the experimental system. 


Actually, the amplitude of the chain oscillation increases 
when the external field is higher and induces larger dipole 
moments in the particles. This increases the tendency of 
the magnetic moments to approach the head-to-tail configu¬ 
ration, which in turn leads to a stronger deformation of the 
matrix. We note that the relative amplitude of the buckling 
along the chain is similar in the simulations (Eig.|^ and exper¬ 
iments (Eig.|^. The matrix surrounding the chain follows the 
chain oscillation with an amplitude that decreases over dis¬ 
tance from the chain. Deviations may be expected from the 
deformational far-held in the experimental system due to the 
different dimensionalities of the systems. 
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FIG. 9: Detailed view of the local deformations in the polymer 
mesh around the magnetic particles with a layer of increased stiff¬ 
ness (top) and without one (bottom) in the immediate vicinity of the 
particle surfaces. The external magnetic field of strength 216 mT is 
applied in the vertical direction. When the boundary layer is assumed 
to be stiffer than the bulk (top), the buckling effect, as observed in the 
experiments, occurs. When the layer around the particles is soft (bot¬ 
tom), neighboring particles either form tight columns parallel to the 
field, or separate in the direction perpendicular to the field. Both, the 
spacing between neighboring particles in the case of a stiff surface 
layer (top) and the overlap of the rigid magnetic particles in the case 
of a soft surface layer (bottom) can be explained by a very strong 
dipole-dipole interaction, which is repulsive in the one case and at¬ 
tractive in the other. 



FIG. 10: Buckling chain of magnetic particles and the surrounding 
polymer mesh for an external field of magnitude 216 mT pointing 
along the vertical direction. In this image, roughly one quarter of the 
full simulation area is shown. The surrounding matrix follows the 
chain oscillation with an amplitude that decreases over distance from 
the chain. 


VI. CONCLUSIONS 

We have shown that paramagnetic chains in a soft polymer 
gel can buckle in a perpendicular magnetic field. The buckling 


morphology depends on the length of the chain, the strength of 
the magnetic field and the modulus of the gel. Longer chains 
form buckling structures with a higher number of half waves. 
Higher strengths of the magnetic field and a lower modulus 
of the gel matrix can lead to higher deformation amplitudes. 
The deformation field in the surrounding gel matrix confirms 
that the embedding polymer network is strongly coupled to 
the paramagnetic chain. A minimal magneto-elastic coupling 
model is developed to describe the morphological behavior of 
the paramagnetic chains in the soft gel under a perpendicular 
magnetic field. It shows that the chains deform in order to de¬ 
crease the magnetic energy. This is hindered by the simultane¬ 
ous deformation of the gel matrix, which increases the elastic 
energy of the gel. Additionally, we have introduced a coarse¬ 
grained molecular dynamics simulation model, which cov¬ 
ers both, the magnetic particles and the surrounding polymer 
mesh. In this model, the buckling of the chains can only be ob¬ 
served when the surface layer around the particles is assumed 
to be stiffer than the bulk of the gel. This prevents the chain 
from breaking up into columns oriented perpendicular to the 
initial chain direction or into isolated particles. These find¬ 
ings support the picture that the embedded magnetic chains 
themselves feature a certain bending rigidity, possibly due to 
polymer chains adsorbed on the particle surfaces. 

Since the magneto-elastic effect demonstrated and analyzed 
in this paper is pronounced, reversible and controllable, it may 
be useful for designing micro-devices, e.g. micro-valves and 
pumps for microfluidic control ||62l. As the morphologies 
of the buckling paramagnetic chains are correlated with the 
modulus of the gel matrix, we may use them as mechanical 
probes for soft gels (similarly to active microrheology tech¬ 
niques) l45l . Moreover, our study may help to understand the 
physical interactions between the magnetic chains and the sur¬ 
rounding cytoskeleton network in magnetotactic bacteria ll46l . 
In our future study we will focus on how the interfacial cou¬ 
pling between the magnetic particles and the polymer network 
influences the local magneto-elastic coupling effect. 
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